Contactless inductive flow tomography 
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The three-dimensional velocity field of a propeller driven liquid metal flow is reconstructed by 
a contactless inductive flow tomography (CIFT). The underlying theory is presented within the 
framework of an integral equation system that governs the magnetic field distribution in a moving 
electrically conducting fluid. For small magnetic Reynolds numbers this integral equation system 
can be cast into a linear inverse problem for the determination of the velocity field from externally 
measured magnetic fields. A robust reconstruction of the large scale velocity field is already achieved 
by applying the external magnetic field alternately in two orthogonal directions and measuring 
the corresponding sets of induced magnetic fields. Kelvin's theorem is exploited to regularize the 
resulting velocity field by using the kinetic energy of the flow as a regularizing functional. The results 
of the new technique are shown to be in satisfactory agreement with ultrasonic measurements. 

PACS numbers: 41.20.Gz, 47.65. +a, 47.80.+V 



I. INTRODUCTION 

Flow measurement in metallic and semiconducting 
melts is a notorious problem in a number of technolo- 
gies, reaching from iron casting to silicon crystal growth. 
Obviously, the usual optical methods of flow measure- 
ment are inappropriate for those opaque fluids. Ultra- 
sonic techniques have problems, too, when applied to 
very hot or chemically aggressive melts. A completely 
contactless flow measurement technique would be highly 
desirable, even if it were only to provide a rough picture 
of the flow. 

Fortunately, metallic and semiconducting melts are 
characterized by a high electrical conductivity. Hence, 
when exposed to an external magnetic field, the flowing 
melt gives rise to electrical currents that lead to a defor- 
mation of the applied magnetic field. This field deforma- 
tion is measurable outside the fluid volume, and it can 
be used to reconstruct the velocity field, quite in paral- 
lel with the well-known magnetoencephalography, where 
neuronal activity in the brain is inferred from magnetic 
field measurements pj . The goal of this paper is to report 
on a first experimental demonstration of such a contact- 
less inductive flow tomography ( CIFT). 



II. THEORY 

The ratio of the induced field to the applied field is 
determined by the so-called magnetic Reynolds number, 
defined as Rm = /iavl, with (i denoting the magnetic 
permeability of the melt, a its electrical conductivity, v 
a typical velocity, and I a typical length scale of the flow. 
In industrial applications, Rm is in the order of 0.01. ..1. 
Only for a few large scale sodium flows, as they appear 
in fast breeder reactors, but also in the recent hydromag- 
netic dynamo experiments |2|, Rm reaches values in the 
order of 10... 100 (of course, in some cosmic dynamos Rm 
can be even much larger) . Actually, the present work was 



strongly motivated by the wish to reconstruct the sodium 
flow in the Riga dynamo experiment by an appropriate 
contactless method. 

Suppose the fluid to flow with the stationary velocity 
v, and to be exposed to a magnetic field B, which we 
leave unspecified for the moment. Then, according to 
Ohms law in moving conductors the current 

j = a(vxB-V^) (1) 

is induced, with ip denoting the electric potential. This 
current gives rise to the induced magnetic field 
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Equation (J2J follows from inserting Eq. (JIJ into Biot- 
Savart's law and transforming the volume integral over 
Vy> into a surface integral over ip. 

The electric potential tp at the boundary S, in turn, 
has to fulfill the boundary integral equation 
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Equation © follows from taking the divergence of Eq. 
(JTJ and utilizing V ■ j = 0. Then, Green's theorem can 
be applied to the solution of the arising Poisson equation 
Aip = V • (v x B), with demanding that the current is 
purely tangential at the boundary 0] ■ Note that Eq. Q 
is the basic formula for the vast area of electric inductive 
flow measurement which is, however, not the subject 
of the present work. 

In general, the magnetic field B on the right hand sides 
of Eqs. 1 ll.il) is the sum of an externally applied mag- 
netic field Bq and the very induced magnetic field b. 
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M - Motor 
P - Propeller 
G - Guiding blades 
H - Helmholtz coils 
B - Hall sensors 




FIG. 1: Scheme (a) and photograph (b) of the CIFT experi- 
ment. 
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FIG. 2: Propeller rotation rate, and induced magnetic field 
measured at the Hall sensor emphasized in Fig. la. 



Hence Eqs. (121: it represent an integral equation system 
which actually can be used to solve dynamo problems in 
arbitrary bounded domains It it also suitable for a 
systematic investigation of the non-linear induction ef- 
fects as they appear already in the sub-critical regime of 
laboratory dynamos [|| . 

In the following, however, all considerations will be 
restricted to problems with small Rm for which B can 
be replaced by Bq. Then, we get a linear relation between 
the desired velocity field and the induced magnetic field 
which is supposed to be measured. But how to cope with 
the remaining Eq. (j3Jl for the electric potential? 

The answer to this question can be adopted from mag- 
netoencephalography Assume, for a given B , all 
measured magnetic field data be collected into an Nb di- 
mensional vector with the entries b\ , and the desired 
velocity components at the Ny discretization points by a 
vector with the entries v n . The solution of the boundary 
integral equation may require a fine discretization of the 
boundary, with Np degrees of freedom (pf°\ Eqs. 100} 
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FIG. 3: Measured induced magnetic field components for 
transverse (a) and axial (b) applied magnetic field, and re- 
constructed velocity as seen from the side (c) and from below 
(d). The grey scale of the arrows indicates the distance from 
the eye. The propeller pumps upward with 1200 rpm. 



can then be written in the form 



R (Bo) v 
t(Bo 
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(4) 
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where the matrices R^ - 1 and T( B °) depend on the ap- 
plied field Bo, whereas the matrices S and U depend only 
on geometric factors. 

As is well known from magnetoencephalography, the 
inversion of Eq. © is a bit tricky due to the singularity of 
the matrix (I — U). This singularity mirrors the fact that 
the electric potential is defined only up to an additive 
constant. We can remove this ambiguity by replacing (I— 
U) by a generally well conditioned matrix (I — TJ) de ^ := 
(I — U) — iV _1 ee T , where e is a vector with all N entries 
equal to one and e T is its transposed. By applying this 
so-called deflation method one ends up with 



'v n + S ik ,(I-U)- k ,Y efl Ttn o) ^ 



(6) 



i.e., with a linear relation between the desired velocity 
field and the measured magnetic field. 

Despite the far-reaching similarity, there is one essen- 
tial difference of our method compared to magnetoen- 
cephalography. While in the latter one has to deter- 
mine a single neuronal current distribution, in our case 
we can produce quite different current distributions from 
the same flow field simply by applying various external 
magnetic fields subsequently. For each applied magnetic 
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FIG. 4: The same as Fig. 3, but for the propeller pumping 
downward with 1200 rpm. 
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FIG. 5: Tikhonov's L-curve for the two different pumping 
directions. The arrows point at the bended knee where the 
curves have the strongest curvature. At these points we get 
a reasonable compromise between data fitting and minimum 
kinetic energy of the modelled velocity field. The rms of these 
velocities is approximately 0.41 m/s for the upward pumping 
and 0.73 m/s for the downward pumping. 



field we can measure the corresponding induced fields, 
and utilize all of them to reconstruct the flow. 

Concerning the uniqueness question for this sort of in- 
version, here we give only a shortened answer, referring 
for more details to the previous papers 0, ■ For spher- 
ical geometry, and the two applied magnetic fields point- 
ing in orthogonal directions, the problem can be solved 
with some rigour. Suppose we have measured the two 



corresponding sets of induced magnetic fields on a sphere 
outside the fluid volume, and have expanded them into 
spherical harmonics. The desired (solenoidal) velocity 
field can be represented by two scalars for its poloidal 
and toroidal parts. These two scalars can also be ex- 
panded into spherical harmonics, but with the expansion 
coefficient still being functions of the radius. In ||| it had 
been shown (at least in some low degrees of the spheri- 
cal harmonics expansion) that what can be derived from 
the two magnetic field expansion coefficients are some ra- 
dial moments of the expansion coefficients for the velocity 
field. A further concretization of the radial dependence 
of the velocity expansion coefficients can only be achieved 
by regularization techniques. If we demand, in a slight 
overinterpretation of Kelvin's theorem, the flow to pos- 
sess minimal kinetic energy, we obtain a unique solution 
for the radial dependence, too. 

Without any rigorous proof at hand, we assume that 
this result can be generalized to aspherical geometry: the 
large main structure of the large scale flow is well infer- 
able, with a depth ambiguity of the velocity that can only 
be resolved by regularization techniques. Imposing two 
orthogonal magnetic fields represents a certain minimum 
configuration for such a flow tomography. For a single 
magnetic field of one direction there are, of course, flow 
components which would be hidden from outside. How- 
ever, all those components are detectable for an external 
magnetic field orthogonal to the previous one. 

For our experimental application we employ the so- 
called Tikhonov regularization [9(, minimizing the total 
functional 



F [v] = F Box [v] + F Boz [v] + F dlv [v] + F reg [v] (7) 
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The first two functionals represent, for applied trans- 
verse field Box and axial field Bo 2 , respectively, the mean 
squared residual deviation of the measured induced mag- 

netic fields b\^ as from the fields b\ [v] modeled ac- 
cording to Eq. ©. Fd^fv] enforces the velocity field to 
be solenoidal, and F reg [v] is the regularization functional 
which tries to minimize the kinetic energy. The param- 
eters (Tj are the assumed a-priori errors for the measure- 
ment of the induced fields. The parameter a c n v is chosen 
very small as it is a measure for the divergence the ve- 
locity solution is allowed to have. The parameter a pen 
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FIG. 6: Axial velocities along the central vertical axis of the 
the cylinder, determined by CIFT and by ultrasonic measure- 
ments (UDV). The ultrasonic measurements are only shown 
up to the propeller position, whereafter they become unreli- 
able. 



determines the trade-off between minimizing the mean 
squared residual deviation of the observed fields and min- 
imizing the kinetic energy of the estimated velocity field. 
The normal equations, that follow from the minimization 
of the functional J7J, are solved by Cholesky decomposi- 
tion. 



III. EXPERIMENT 

In the experiment (Fig. 1) we use 4.4 liters of the eutec- 
tic alloy Ga 67 In 20 5 Sn 12 - 5 that is liquid at room temper- 
atures. The flow is produced by a motor driven propeller 
with a diameter of 6 cm inside a cylindrical polypropylene 
vessel with 18.0 cm diameter. The height of the liquid 
metal is 17.2 cm, yielding an aspect ratio close to 1. 

The position of the propeller is approximately at one 
third of the total hight, measured from the top. Eight 
guiding blades above the propeller are intended to re- 
move the swirl of the flow for the case that the propeller 
pumps upward. Contrary to that, the downward pump- 
ing produces, in addition to the main poloidal motion, 
a considerable toroidal motion. The rotation rate of the 
propeller can reach up to 2000 rpm, which amounts to a 
mean velocity of approximately 1 m/s, corresponding to 
a magnetic Reynolds number of approximately 0.4. 

Two pairs of Helmholtz coils are fed by currents of 22.5 
Ampere and 32.5 Ampere, respectively, to produce alter- 
nately an axial and a transversal field of 4 mT, which 
both are rather homogeneous throughout the vessel. Ei- 
ther field is applied for a period of 3 seconds, during 
which a trapezoidal signal form is used. The measure- 
ments are carried out for 0.5 seconds, 1 second after the 
plateau value of the trapezoidal current has been reached. 
Hence, we get an online monitoring with a time resolu- 
tion of 6 seconds. 




FIG. 7: Sketch of the ultrasonic measurement set-up for the 
azimuthal velocity component. The obtained velocity along 
the chord is a position dependent mixture of radial and az- 
imuthal components. In the middle of the chord, one gets the 
pure azimuthal component, (a) Typical situation for down- 
ward pumping with dominant azimuthal velocity (at the axial 
position 70 mm from the bottom). The projection of the ve- 
locity onto the chord has a maximum in the middle of the 
chord, (b) Typical situation for upward pumping with dom- 
inant radial velocity. The projection of the velocity goes to 
zero in the middle of the chord. 



The induced magnetic fields are measured by 49 Hall 
sensors, 8 of them grouped together on each of 6 cir- 
cuit boards which are located on different heights (Fig. 
1). One additional sensor is located in the center below 
the vessel. The key problem of the method is the re- 
liable determination of comparably small induced mag- 
netic fields on the background of much higher imposed 
magnetic fields. An accurate control of the external mag- 
netic field is essential to meet this goal. In our configura- 
tion the current drift in the Helmholtz coils can be con- 
trolled with an accuracy of better than 0.1 per cent. This 
is sufficient since the measured induced fields are approx- 
imately 1 per cent of the applied field. The temperature 
drift of the sensitivity can be overcome by enforcing the 
applied current in the Hall sensors to be constant. The 
temperature drift of the offset problem is circumvented 
by changing the sign of the applied magnetic field. Fig- 
ure 2 shows that by these means a stable measurement 
of the small induced field can be realized, even over a 
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FIG. 8: Measured projection of the velocity onto the chord, 
for downward and upward pumping. In the middle of the 
chord (at 85 mm distance from the wall) we get an azimuthal 
velocity of +0.58 m/s for downward pumping and -0.05 m/s 
for upward pumping (note the change of sign due to con- 
ventions in UDV measurements). Compare also Fig. 7 for 
illustration. 
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FIG. 9: Azimuthal velocity at r=30 mm, as determined by 
CIFT at different axial positions. The small symbols (crosses 
and squares) represent individual measurements at six az- 
imuthal positions. The large crosses and squares represent 
the corresponding averages over six azimuthal positions. The 
full symbols (circle and triangle) are the UDV values as in- 
ferred from Fig. 8. 



period of one hour. 

For upward and downward pumping, Figs. 3 and 4 
show the induced magnetic fields measured at the 49 po- 
sitions, and the inferred velocity field at 52 discretization 
points. In Fig. 3c we see clearly the upward flow in the 
center of the vessel and the downward flow at the rim, 
but nearly no rotation of the flow in Fig. 3d. In Fig. 4c 
we can identify the downward flow in the center and the 
upward flow at the rim, and in Fig. 4d a clear rotation 
of the flow. Evidently, the method is able to identify the 



poloidal rolls and the absence or presence of the swirl. 

For both flow directions, Fig. 5 illustrates the applica- 
tion of Tikhonov's L-curves [9j. This curve, which results 
from scaling the parameter a pen in Eq. (11) from lower to 
higher values, shows the dependence of the mean squared 
residual of the measured data on the kinetic energy of the 
flow. For low values (left end of Fig. 5) only little kinetic 
energy is allowed, leading to a velocity field that fits the 
measured magnetic field data only poorly. For high val- 
ues (right end of Fig. 5) the data are fitted very well but 
with an unphysical high kinetic energy. At the points of 
strongest bending (the "knee"), the resulting velocities 
(Figs. 3 and 4) are physically most reasonable 9]. 



IV. VALIDATION 

In order to validate the CIFT method, we have per- 
formed independent velocity measurements based on ul- 
trasonic Doppler velocimetry (UDV). For that purpose 
we have used the DOP2000 ultrasonic velocimeter man- 
ufactured by Signal-Processing SA (Lausanne, Switzer- 
land), which had already demonstrated its capabilities 
for velocity measurements in liquid metals |lfl Hlj . As 
ultrasonic transducers we have used 2 MHz probes. 

Because of its comparably large magnetic Reynolds 
number (Rm « 0.2), the propeller driven flow in the 
cylinder has also a large hydrodynamic Reynolds num- 
ber (Re w2x 10 5 ). Necessarily, the flow is highly turbu- 
lent. Strong fluctuations are observed both by the CIFT 
method as well as by UDV. 

For a sensible comparison of both methods, some time 
averaging is advised. In the following we will focus on two 
UDV measurements that were both taken at a propeller 
rotation rate of 1200 rpm, and which represent a time 
average over half a minute. 

The first measurement concerns the axial velocity 
along the central vertical axis of the cylinder. This axial 
component is easily measured by an ultrasonic transducer 
flash mounted to the bottom of the cylinder. Figure 6 
shows the results of the UDV measurement (up to the 
propeller position), together with the results of the CIFT 
measurement. For both upward and downward pump- 
ing we see a reasonable correspondence of both measure- 
ments. Notably, CIFT exhibits the different axial de- 
pendencies that are typical for upward and downward 
pumping and which are confirmed by the UDV data. 

The second measurement, which concerns the az- 
imuthal velocity component, deserves some explanation. 
Figure 7 shows the UDV measurement set-up. The axial 
position is at 70 mm from the bottom. What is actu- 
ally measured by UDV is the projection of the velocity 
onto the ultrasound beam along the chord. Therefore, 
the measured signal is in general a mixture of the radial 
and azimuthal velocity components. Only in the middle 
of the chord we get a signal that originates purely from 
the azimuthal velocity In Figs. 7 (a) and (b) we illus- 
trate the measured data that are shown in Fig. 8. In 
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the case of downward pumping the velocity is dominated 
by the rotation whereas for upward pumping it is domi- 
nated by the radial part. In the middle of the chord we 
infer a mean azimuthal velocity of 0.58 m/s for downward 
pumping and -0.05 m/s for upward pumping. 

Do these UDV values agree with those obtained by 
CIFT? In Fig. 9 we show the axial dependence of the 
azimuthal velocity at a radial position r=30 mm, as ob- 
tained by CIFT during a measurement time of 6 sec- 
onds. We give here the individual values at six different 
azimuthal positions (small crosses and squares). Inter- 
estingly, though the system is in general axisymmctric, 
non-axisymmetric fluctuations are still visible. The av- 
erages (large crosses and squares) over the six azimuthal 
positions show, at an axial position of z=70 mm, a good 
agreement with the data from UDV measurement. 

V. CONCLUSIONS AND PROSPECTS 

To summarize, we have put into practice a first version 
of contactless inductive flow tomography, using two or- 
thogonal imposed magnetic fields. The comparison with 
UDV measurements shows that the method provides ro- 
bust results on the main structure and the amplitude 
of the velocity field. A particular power of CIFT con- 
sists in a transient resolution of the full three-dimensional 
flow structure in steps of several seconds. Hence, slowly 



changing flow fields in various processes can be followed 
in time. Due to its weakness the externally applied mag- 
netic field does not influence the flow to be measured. 
However, CIFT is also possible in cases where stronger 
magnetic fields are already present for the purpose of flow 
control, as, e.g., the electromagnetic brake in steel cast- 
ing or the DC-field components in silicon crystal growth. 
Obviously, the future of the method lays with applying 
AC fields with different frequencies in order to improve 
the depth resolution of the velocity field. For problems 
with higher Rm, including dynamos, the inverse prob- 
lem becomes non-linear, and more sophisticated inversion 
methods must be applied to infer the velocity structure 
from magnetic field data. Although interesting results 
have been obtained by employing Evolutionary Strate- 
gies to inverse spectral dynamo problems [12|, and first 
tests of such inversion schemes for the data from the Riga 
dynamo experiment have shown promising results, the 
general inverse dynamo topic is extremely complicated 
and goes essentially beyond the scope of the present pa- 
per. 
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